Condensate formation with three-component ultracold fermions 
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We investigate the formation of Bose-Einstein condensation and population imbalance in a three- 
component Fermi superfluid by increasing the U(3) invariant attractive interaction. We consider 
the system at zero temperature in three dimensions and also in two dimensions. Within the mean- 
field theory, we derive explicit formulas for number densities, gap order parameter, condensate 
density and condensate fraction of the uniform system, and analyze them in the crossover from the 
Bardeen-Cooper-Schrieffer (BCS) state of Cooper pairs to the Bose-Einstein Condensate (BEC) of 
strongly-bound molecular dimers. In addition, we study this Fermi mixture trapped by a harmonic 
potential: we calculate the density profiles of the three components and the condensate density 
profile of Cooper pairs in the BCS-BEC crossover. 
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I. INTRODUCTION 

In the last years two experimental groups P, [1] 
have analyzed the condensate fraction of ultra-cold two 
hyperfine component Fermi vapors of ^Li atoms in 
the crossover from the Bardeen-Cooper-SchriefFer (BCS) 
state of Cooper Fermi pairs to the Bose-Einstein conden- 
sate (BEC) of molecular dimers. These experiments are 
in quite good agreement with mean-field theoretical pre- 
dictions [3, 4] and Monte-Carlo simulations 5] at zero 
temperature, while at finite temperature beyond-mean- 
field corrections are needed (Q] . 

Degenerate three-component gases have been experi- 
mentally realized using the three lowest hyperfine states 
of ^Li 0, Q. At high magnetic fields the scattering 
lengths of this three-component system are very close 
each other and the system is approximately U(3) in- 
variant. Moreover, it has been theoretically predicted 
that good SU{N) invariance (with N < 10) can be 
reached with ultracold alkaline-earth atoms (e.g. with 
^''Sr atoms) Very recently, Ozawa and Baym have 
investigated the BCS-BEC crossover of uniform three- 
component ultracold fermions with U(3) symmetry at 
zero and finite temperature in three-dimensional space 

In this paper we focus on the condensate formation 
in a uniform three-component Fermi superfluid by in- 
creasing the U(3) invariant attractive interaction. We 
determine the condensate fraction and the population 
imbalance of the system in both the three-dimensional 
case and in the two-dimensional one. By using the ex- 
tended BEC theory [ll|, we obtain analytical formulas 
for number densities, condensate density and condensate 
fraction 0,[l3 of the uniform system and analyze them in 
the crossover from the Bardeen-Cooper-SchriefFer (BCS) 
state of Cooper pairs to the Bose-Einstein Condensate 
(BEC) of strongly-bound molecular dimers. We consider 
also the inclusion of a harmonic confinement, determining 
the density profiles of the three components of the ultra- 
cold gas, and the space dependence of condensate density 



and gap order parameter. In the paper we neglect effects 
of three-body interactions like Efimov three-body bound 
states. 



II. THE MODEL 

The shifted Hamiltonian density of a dilute and inter- 
acting three-component Fermi gas in a volume V is given 
by [10] 



2m 



a=R,G,B 
+ i^R'^'ti^B^R + 'IpG^t'ipB'IpG) 



(1) 



where ipa (r) is the field operator that destroys a fermion 
of component a in the position r, while i/'^(r) creates a 
fermion of component a in r. To mimic QCD the three 
components are thought as three colors: red (R), green 
(G) and blue (B). The attractive inter-atomic interaction 
is described by a contact pseudo-potential of strength g 
{g < 0). The number density operator is 



n(r) = ^ V)+(r)?/>c,(r) 

a=R,G,B 



and the average number of fermions reads 



N = 



r (n 



(r)). 



(2) 



(3) 



This total number N is fixed by the chemical potential 
/J, which appears in Eq. ([TJ. As stressed in Ref. [lo| . 
by fixing only the total chemical potential ^ (or equiva- 
lently only the total number of atoms N) the Hamilto- 
nian ([1]) is invariant under global U(3) rotations of the 
species (and consequently also under global SU(3)). At 
zero temperature, the attractive interaction leads to pair- 
ing of fermions which breaks the U(3) symmetry but only 
two colors are paired and one is left unpaired [lOj. 
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Following Ozawa and Baym we assume without 
loss of generality that the red and green particles are 
paired and the blue are not paired. The interacting terms 
can be then treated within the minimal mean-field BCS 
approximation, giving 

g ■4;^'ip^'ipc'ipR = g {'ip%i)Q)^Gi^R + 9 ^R^t:{'iGi'R) (4) 



and 



^R^B^BlpR = 9 i)%i)'^i)Bi^G = . 



(5) 



Notice that the Hartree terms have been neglected, while 
the pairing gap A — g{iipG'>pR} between red and green 
fcrmions is the key quantity. 

The shifted Hamiltonian density ([IJ is diagonalized by 
using the Bogoliubov-Valatin representation of the field 
operator ipaii^) in terms of the anticommuting quasi- 
particle Bogoliubov operators bka with amplitudes Uk 
and Vk and energy Ek ■ After minimization of the result- 
ing quadratic Hamiltonian one finds familiar expressions 
for these quantities: 
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= ^ 1 - 
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Ek 



(6) 



(7) 



with u\ — 1 — v^.. In addition wc find the equation for 
the number of particles 



N = Nr + Ng + Nb 



where 



and 



Nr = Ng^\Y.''^ 



NB^Y^^iy^- 
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(8) 



(9) 



(10) 



with Q{x) the Heaviside step function, and also the gap 
equation 



1 _ 1 ^ 1 
5 ^V^2E^ 



(11) 



The chemical potential fi and the gap energy A are ob- 
tained by solving equations (|8]) and pT|) . 

We observe that the number of red-green pairs in the 
lowest state, i.e. the condensate number of red-green 
pairs, is given by d, |3, Q 

TVo = / dh, |(V'G(ri)^i?(r2))P, (12) 



and it is straightforward to show that 



^^0 = $^M^ 



2 

k ■ 



(13) 



In the continuum limit, due to the choice of a contact 
potential, the gap equation (jlip diverges in the ultravi- 
olet. This divergence is linear in three dimensions and 
logarithmic in two dimensions. Let us face this problem 
in the next two sections. 



III. THREE-COMPONENT ULTRACOLD 
FERMIONS IN THREE DIMENSIONS 

In three dimensions, a suitable regularization [11] is ob- 
tained by introducing the inter-atomic scattering length 
ap via the equation 
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E 



m 



(14) 



and then subtracting this equation from the gap equa- 
tion (jlip . In this way one obtains the three-dimensional 
regularized gap equation 
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V^\2Ek h'^k'^ 



(15) 



In the three-dimensional continuum limit — >■ 
V/{2'kY J d^h y/(27r^) / k'^dk from the number equa- 
tion ([SJ with (jni) and (ITUl) we find the total number den- 
sity as 



n = y = riR + nc + riB : 



with 



and 



1 (2m)3/2 / u 



(16) 



(17) 
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3 27T^h^ 

The renormalized gap equation (I15|) becomes instead 



(18) 



1 ^ 2(2m)V2 /^x 



(19) 



where kp = {<oi:N / {W)Y/^ = {2Tr'^n)^/^ is the Fermi 
wave number. Here Ii{x) and l2ix) are the two mono- 
tonic functions 



hix) = / y 



+ 00 



1 



1 




dy , (20) 



dy , (21) 
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which can be expressed in terms of elhptic integrals, as 
shown by Marini, Pistolesi and Strinati In a similar 
way we get the condensate density of the red-green pair 

as 



no = 



No 
V 



^3/2 




(22) 



This equation and the gap equation ()19|) are the same 
of the two-component superfluid fermi gas (see 0]) but 
the number equation ([8]), with (|9]) and ([TOl). is clearly 
different. Note that all the relevant quantities can be 
expressed in terms of the ratio 



Xq 



A 



(23) 



For instance, the fraction of red fermions, which is equal 
to the fraction of green fermions, is given by 

rvR _nG_ _ h{xn) 

n n 



2h{xo) + Ixl'^ e(a;o) 



(24) 



while the fraction of blue fermions reads 



1 



n 



The fraction of condensed red-green pairs is instead 



(25) 



no 
n 



\Jxo + y/l + xl 



h{xo) + \xl'^ e(xo) 



(26) 



Finally, the adimensional interaction strength of the 
BCS-BEC crossover is given by 

1 2 h{xo) 



y 



kpap 



^2h{xo) 



2^3/2 p, 

3 ' 



Q{xo) 



(27) 



We can use these parametric formulas of to plot the 
density fractions as a function of the scaled interaction 
strength y. 

In the upper panel of Fig. [T] we plot the fraction of 
red fermions n^ jn (solid line) and the fraction of blue 
fermions UB/n (dashed line) as a function of scaled in- 
teraction strength y = 1/ (kpap)- The behavior of UQ/n 
is not shown because it is exactly the same of nn/n. The 
figure shows that in the deep BCS regime (y <C —1) the 
system has npi/n — ndn = UB/n = 1/3. By increas- 
ing y the fraction of red and green fermions increases 
while the fraction of blue fermions decreases. At j/ ~ 0.6, 
where xq — 0, the fraction of blue fermions becomes zero, 
i.e. ub/ti = and consequently nji/n — uq/ii — 1/2. 
For larger values of y there are only the paired red and 
green particles. This behavior is fully consistent with the 
findings of Ozawa and Baym [l^ . 

In the lower panel of Fig (T] it is shown the plot of the 
condensate fraction no/{n/2) of red-green pairs through 
the BCS-BEC crossover as a function of the Fermi-gas 
parameter y = l/ikpaF)- The figure shows that a large 
condensate fraction builds up in the BCS side already 
before the unitarity limit (j/ = 0), and that on the BEC 
side (y ^ 1) it rapidly converges to one. 
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FIG. 1: (Color online). Three-component ultracold fermions 
in three-dimensions. Upper panel: fraction of red fermions 
jifl/n (solid line) and fraction of blue fermions ns/n (dashed 
line) as a function of scaled interaction strength y = 
l/ikpaF)- Lower panel: condensed fraction of red-green 
particles no/n as a function of scaled interaction strength 
y = l/{kFaF). 



IV. THREE-COMPONENT ULTRACOLD 
FERMIONS IN TWO DIMENSIONS 

A two-dimensional Fermi gas can be obtained by im- 
posing a very strong confinement along one of the three 
spatial directions. In practice, the potential energy Ep 
of this strong external confinement must be much larger 
than the total chemical potential h^d of the fermionic 
system: /is^i ^ 2Ep [ij. Contrary to the three- 
dimensional case, in two dimensions quite generally a 
bound-state energy cb exists for any value of the inter- 
action strength g between atoms [li|. For the contact 
potential the bound-state equation is 



1 



1 



^ k 2m ^ 



(28) 



and then subtracting this equation from the gap equa- 
tion (llip one obtains the two-dimensional regularized gap 
equation [ll| 



k 



2m 



= 



(29) 



Note that, for a 2D inter-atomic potential described 
by a 2D circularly symmetric well of radius i?o and 
depth Uo, the bound-state energy es is given hy tB — 
h^/{2mRl)eyi]y{^2h^/{mUoRl)) with U^Rl ^ [H. 
In the two-dimensional continuum limit ^j^. — )■ 
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FIG. 2: (Color online). Three-component ultracold fermions 
in two-dimensions. Upper panel: fraction of red fermions 
HR/n (solid line) and fraction of blue fermions UB/n (dashed 
line) as a function of scaled bound-state energy esj^F ■ Lower 
panel: condensed fraction of red-green particles no/n as a 
function of scaled bound-state energy tsl^F- 



V/{2'!Tf J d^k V/{27r) J kdk, the Eq. ^ gives 




(30) 



Instead, the number equation ^ with ^ and ([TU)) gives 
the total number density as 



N 
V 



nR + nG + ub 



(31) 



where is a two-dimensional volume (i.e. an area), the 
red and green densities are 



while the blue density reads 




(32) 



(33) 



Finally, the condensate density of red-green pairs is 
given by 



1 , m 



no 



4 ^ 271^2 



)A(|-f arctan(£)) 



(34) 



Also in this two-dimensional case all the relevant quanti- 
ties can be expressed in terms of the ratio xq = fJ./A. In 



particular, the fraction of red fermions, which is equal to 
the fraction of green fermions, is given by 



riR 
n 



n 



2[xo + ^l + xl + XQ Q{xo)] 



(35) 



the fraction of blue fermions is 



- 1 _ 2 — 



and the condensate fraction is 



(36) 



no 
n 



^ -I- arctan (xq) 



4[xo + yr+ai + 2:0 e(a;o)] 



(37) 



In two-dimensions the Fermi energy ep = h'^kp/{2m,) of 
a non-interacting Fermi gas is given by e^? = irfi^n/m 
with kp = iinN/{3V)y/^ = (47rn/3)i/2 the Fermi wave 
number. It is convenient to express the bound-state en- 
ergy £5 in terms of the Fermi energy ep. In this way we 
find 



y/l + xl- Xo 



Xq + ^/T+xl + Xo 6(a;o) 



(38) 



We can now use these parametric formulas of xq to plot 
the fractions as a function of the scaled bound-state en- 
ergy es/ep- 

In the upper panel of Fig. [5] we plot the fraction 
of red fermions nn/n (solid line) and the fraction of 
blue fermions ns/n (dashed line) as a function of scaled 
bound-state energy es/^F- The behavior of uqIji is not 
shown because it is exactly the same of nj^/n. The fig- 
ure shows that in the deep BCS regime {eB/(-F ^ 1) 
the system has n^ jn — UQ/n = Nb/tl = 1/3. By in- 
creasing cb/ the fraction of red and green fermions in- 
creases while the fraction of blue fermions decreases. At 
^b/^f = 3, where = 0, the fraction of blue fermions 
becomes zero. For larger values of €b/^f there are only 
the paired red and green particles. This behavior is quite 
similar to the one of the three-dimensional case; the main 
difference is due to the fact that here the curves are lin- 
ear. 

In the lower panel of FigH] it is shown the conden- 
sate fraction tiq/ {n/2) of red-green pairs. In the weakly- 
bound BCS regime {tB/^F ^ 1) the condensed frac- 
tion no/n goes to zero, while in the strongly-bound BEC 
regime (^bI^f ^ 1) the condensed fraction no/^i goes 
to 1/2, i.e. all the red-green Fermi pairs belong to the 
Bose-Einstein condensate. Notice that the condensate 
fraction is zero when the bound-state energy es is zero. 
For small values of es/f-p the condensed fraction has a 
very fast grow but then it reaches the asymptotic value 
1/2 very slowly. 
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FIG. 3: (Color online). Three-component ultracold fermions 
in two-dimensions under harmonic confinement U{r) — 
Left panels: density profile n_R(r) of red fermions (solid lines) 
and density profile nsir) of blue fermions (dashed lines). 
Right panels: total density profile n(r) = 2nR{r) + nB{r) 
(solid lines) and condensate density profile 2no (r) (dot-dashed 
lines). Note that the density ncir) of green fermions is not 
shown because it is equal to nair). Data obtained with 
Jj, = IQ and three values of the bound-state energy es: e_B = 1 
(upper panels), es = 10 (middle panels), and es = 50 (lower 
panels). We use scaled variables where lengths are in units of 
an = y/KJJrmjj), with uj the frequency of harmonic confine- 
ment, and densities in units of aT? . 



INCLUSION OF HARMONIC 
CONFINEMENT 



blue fermions read 



nB{r) = i^)Kr) ^iKr)) 




(43) 



The density profile of condensed red-geen pairs is instead 
given by 



arctan 



A(r)- 



(44) 



In Fig. [3] we plot results obtained by using scaled 
variables; energy in units of fku and length in units of 
an = \/h/ {niLo). We work at fixed chemical potential 
/i = 10 and increase the bound-state energy e^- In the 
upper panels of Fig. |3] we set €b = 1, in the middle 
panels es = 10 and in the lower ones = 50. In the 
left panels of Fig. [3] we show density profiles of red par- 
ticles (solid lines) and blue particles (dashed lines). In 
the right panels we show instead the total density pro- 
file n(r) = ni^{r) -\- ncir) + nsir) (solid line) and the 
density profile 2no(r) of condensed particles (dot-dashed 
lines). As expected, by increasing ts the critical radius 
To = y/^B + 2/i (written in scaled units), at which the gap 
order parameter A(r) and the condensate density no(r') 
go to zero, becomes much larger than the critical radius 
tt = -^72^, at which the density Ub (r) of blue fermions 
goes to zero. Moreover, by increasing cb the cloud of 
blue particles becomes quite small while the cloud of con- 
densed particles approaches the total one. 



It is interesting to study the effect of a harmonic po- 
tential 

U{r) = ^muj^r"^ (39) 

on the properties of the three-component ultracold gas 
in the BCS-BEC crossover. For semplicity we investigate 
the two-dimensional case, which gives rise to elegant for- 
mulas also in this non-uniform configurtion. In fact, by 
using the local density approximation, namely the sub- 
stitution 

^l ^ ^i{r) - p - U[r) , (40) 

the gap equation (pO| gives the space-dependent gap pa- 
rameter as 

2 2 

A(r) = Ao (1-^)6(1-^), (41) 

where Aq — y^e^ -|- 2eB/i and ro = Aq/ ^/e^rnUj. Here p, 
is the chemical potential of the non-uniform system. 
In the same way the density profiles of red, green and 



VI. CONCLUSIONS 

We have investigated the condensate fraction and the 
population imbalance of a uniform three-component ul- 
tracold fermions by increasing the U(3) invariant attrac- 
tive interaction. We have considered the superfiuid sys- 
tem both in the three-dimensional case and in the two- 
dimensional one. For the uniform system we have ob- 
tained explicit formulas and plots for number densities, 
condensate density and population imbalance in the full 
BCS-BEC crossover. For the system under harmonic 
confinement we have analyzed the density profiles of the 
three components and the density profile of the con- 
densed pairs by varying the interaction strength. In our 
calculations we have used the mean-field extended BCS 
theory. Monte Carlo simulations have shown that, at 
zero-temperature, beyond-mean-field effects are negligi- 
ble in the BCS side of the BCS-BEC crossover while they 
become relevant in the BEC side [1, [ll]. Nevertheless, in 
the deep BEC side the condensate fraction goes in any 
case to one and the main difference in its determination 
is around the unitarity limit, where the relative difference 
between mean-field and Monte Carlo condensate fraction 



6 



is about 20% [H . We think that our results can be of in- 
terest for next future experiments with degenerate gases 
made of alkah-metal or alkaline-earth atoms. As stressed 
in the introduction, SU(N) invariant interactions can be 
experimentaly obtained by using these atomic species 0- 
1^. The problem of unequal couplings, and also that of a 



fixed number of atoms for each component, is clearly of 
big interest too, and its analysis can be afforded by in- 
cluding more than one order parameter (see for instance 

M)- 

The author thanks Flavio Toigo for useful discussions 
and Nicola Manini for computational help. 
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